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We study the influence of cooperative effects such as superradiance and subradiance, on the 
scattering properties of dilute atomic gases. We show that cooperative effects lead to an effective 
potential between pairs of atoms that decays like 1/r. In the case of superradiance, this potential 
is attractive for close enough atoms and can be interpreted as a coherent mesoscopic effect. We 
consider a model of multiple scattering of a photon among superradiant pairs and calculate the 
elastic mean free path and the group velocity. We study first the case of a scalar wave which allows 
to obtain and to understand basic features of cooperative effects and multiple scattering. We then 
turn to the general problem of a vector wave. In both cases, we obtain qualitatively similar results 
and derive, for the case of a scalar wave, analytic expressions of the elastic mean free path and of 
the group velocity for an arbitrary (near resonance) detuning. 
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I. INTRODUCTION 

Coherent multiple scattering of photons in cold atomic 
gases is an important problem since it presents a path 
towards the onset of Anderson localization transition, a 
long standing and still open issue. The large resonant 
scattering cross-section of photons reduces the elastic 
mean free path to values comparable to the photon wave- 
length for which the weak disorder approximation breaks 
down thus signaling the onset of Anderson localization 
transition [1, 2]. Another advantage of cold atomic gases 
is that sources of decoherencc and inelastic scattering 
such as Doppler broadening are often negligible. More- 
over, propagation of photons in atomic gases differs from 
the case of electrons in disordered metals or of electro- 
magnetic waves in suspensions of classical scatterers for 
which mesoscopic effects and Anderson localization have 
been thoroughly investigated [1]. This problem is then 
of great interest since it may raise new issues in the An- 
derson problem such as change of universality class and 
therefore new critical behavior. New features displayed 
by the photon-atom problem are the existence of internal 
degrees of freedom (Zecman sublevels) and cooperative 
effects such as subradiance or superradiance [3] that lead 
to effective interactions between atoms [4]. These two 
differences are expected to lead to qualitative changes 
of both mesoscopic quantities and Anderson localization. 
Some of the effects of a Zeeman degeneracy have been 
investigated in the weak disorder limit [5] using a set of 
finite phase coherence times [6] which reduce mesoscopic 
effects, such as coherent backscattering [1, 7]. 

The influence of cooperative effects and more specifi- 
cally of superradiance on the multiple scattering of pho- 
tons has been recently investigated [8] . It has been shown 
that in atomic gases superradiance and subradiance lead 
to a potential between two atoms, analogous to the one 
considered in [9, 10], that decays like the inverse of the 
distance between them. In the case of superradiance, 



this potential is attractive for close enough atoms, and 
can be interpreted as a coherent mesoscopic effect. The 
contribution of superradiant pairs to multiple scattering 
properties of a dilute gas has been calculated by using 
an effective propagator that describes a scalar wave being 
scattered by a pair of two-level atoms. Simple expressions 
for the photon elastic mean free path and group velocity 
have been derived at resonance and found to be signifi- 
cantly different from that of independent atoms. To be 
more specific, near resonance, as well as at resonance, the 
superradiant effect leads to a finite and positive group ve- 
locity, unlike the one obtained for light interaction with 
independent atoms. 

In this paper we provide, for the case of a scalar wave, 
closed expressions for the suprerradiant contribution to 
the elastic mean free path and the group velocity for an 
arbitrary (near resonance) detuning, and calculate the 
dependence of the transport time on it. In addition, we 
estimate the maximal inter-atomic separation in a su- 
perradiant pair, which accounts for possible mechanisms 
that may break the pair. We also compare the effective 
approach presented in [8] to a more realistic one that 
takes into account the vectorial nature of the wave. 

The paper is organized as follows: In section II we 
describe the model which consists of pairs of two-level 
atoms placed in an external radiation field where the 
Doppler shift and recoil effects are negligible. In order 
to investigate the influence of the cooperative effects of 
such pairs on the multiple scattering of photons we briefly 
review, in section III, Dickc states and some of their prop- 
erties. Then, we calculate the average interaction poten- 
tial of a pair of atoms in a Dicke state by averaging upon 
the random orientations of pairs of atoms with respect to 
the wave vector of a photon incident on the atomic cloud. 
Next, we study the scattering of a photon by such pairs 
and, in section IV, compare the results to the case where 
a classical wave is being scattered by a pair of atoms. 
This comparison allows to find an unexpected connec- 
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tion between superradiance and mcsoscopic effects. In 
sections V and VI, we consider the multiple scattering of 
photons by pairs of atoms and calculate the elastic mean 
free path and the group velocity of photons in the ran- 
dom medium. Finally our analysis is compared to other 
approaches in part VII and its results are summarized in 
part VIII. 

II. MODEL 



We assume that the typical speed of the atoms, v ~ 
■^/fcsTo/jU, is small compared to v rnax ~ T/k but large 
compared to v m i n — hkj ' [i where fi is the mass of the 
atom and To is the temperature, so that it is possible 
to neglect the Doppler shift and recoil effects. Indeed, 
for a temperature of To = lO -3 ^, the typical speed of 
the atom is v ~ 0.3m/ s. Since, for a wave number of 
k = 10 7 m _1 and a natural width of T = 10 7 s _1 , v max ~ 
lm/ s and v m i n ~ O.Olm/s both assumptions are fulfilled. 



Atoms are taken to be degenerate, two-level systems 
denoted by \g) = \j g = 0, m g = 0) for the ground state 
and |e) = |j e = l,m e = 0, ±1) for the excited state, 
where j is the total angular momentum and m is its pro- 
jection on a quantization axis, taken as the z axis. The 
energy separation between the two levels, including ra- 
diative shift, is fuvo and the natural width of the excited 
level is KT. This simple picture of a two- level atom ne- 
glects the rather complicated energy structure of a real 
atom which reflects various internal interactions, e.g., 
Coulomb interactions, spin-orbit interactions, hypcrfmc 
interactions, etc. But, due to selection rules which limit 
the allowed transitions between states, in some cases a 
certain state may couple to only one other. Thus, the 
two-level atom approximation is close to reality and not 
merely a mathematical convenience. 

We consider a pair of such atoms in an external ra- 
diation field and the corresponding Hamiltonian is H = 
H a + V, with 

h» = ^£>)<e| - \B)(9\)i + £^ a L«k £ . (1) 

1=1 ke 

ak e (resp. ) is the annihilation (resp. creation) oper- 
ator of a mode of the field of wave vector k, polarization 
£k and angular frequency cuk = c|k|. The interaction V 
between the radiation field and the dipole moments of 
the atoms is given by 

2 

V = -^d ; .E(r ; ) (2) 

i=i 

where E(r) is the electric field operator 

E(r) = i J2 {^^^' T ~ 4e4e- ik ' r ). (3) 

ke 

f2 is a quantization volume and d; is the electric dipole 
moment operator of the Z-th atom. As an odd operator, 
which changes sign upon inversion, d; may be written as 

d ; = ( 5 |d|e)Ar + (e|d| fl )A+ (4) 

where the atomic raising and lowering operators are 

A+ = (\e)(g\) l A," = (IffXel), . (5) 



III. DICKE STATES 
A. Interaction potential and lifetime 

The absorption of a photon by a pair of atoms in 
their ground state leads to a configuration where the two 
atoms, one excited and the second in its ground state, 
have multiple exchange of a photon, giving rise to an ef- 
fective interaction potential and to a modified lifetime as 
compared to independent atoms. These two quantities 
are obtained from the matrix elements of the evolution 
operator U(t) between states such as |gie 2 ;0). There 
are six unperturbed and degenerate states with no pho- 
ton, given by {|<7ie2i; 0), |eij<72; 0}} in a standard basis 
where i,j = —1,0,1. The symmetries of the Hamilto- 
nian, namely its invariance by rotation around the axis 
between the two atoms, and by reflection with respect to 
a plane containing this axis, allows to use combinations 
of these states that are given by 

\4>!) = ^=[|ei i52 ;0}+e| 5l e 2i ;0)] (6) 

with e = ±1, so that 

(d><'\U(tM) = 6 ij 5 ee ,Sm (7) 

and 

Sl(t) = (e ll g2;0\U{t)\e ll g 2 -,0)+e{g 1 e 2l -,0\U{t)\e 1 , l g 2] 0). 

(8) 

The states \<j>f) may be rewritten in terms of the well- 
known Dicke states \LM), where L is the cooperation 
number and M is half of the total atomic inversion [3]. 
For two atoms, the singlet Dicke state is 

|00) = -^[\em) - \91e2)] (9) 

and the triplet Dicke states are 
|H> = |eie 2 ) 
|10> = ^=[|ei<fc> + |ffie 2 >] 

|1-1) = \9i9%) (10) 

The states |11) and |1 — 1) correspond, respectively, to 
both atoms in their excited states and both atoms in 
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their ground state. The singlet state |00) and the triplet 
state 1 10) both correspond to one atom in the excited 
state and the other in the ground state, but |00) is anti- 
symmetric where |10) is symmetric under an exchange 
of the atoms. Therefore, we may rewrite (6) as \<f>f) = 
|10;0) and \4>r) = |00;0). 

For times such that t 3> r/c, where r is the distance be- 
tween the two atoms, up to second order in the coupling 
to the radiation, (8) reads 



i 



.fin 



fin 



The two real quantities AEf and r| are respectively the 
interacting potential and the probability per unit time of 
emission of a photon by the two atoms in the state \4>\). 
The calculation of these two quantities requires second 
order perturbation theory with respect to the interaction 
(2). To this purpose we define an initial state where 
one atom is excited and the other is in its ground state 
without any photon and a final state where the two atoms 
are exchanged. We also define intermediate states of two 
types: both atoms in their ground state with one virtual 
photon present and both atoms in their excited state with 
one virtual photon present. Summing the corresponding 
diagrams [11] gives 



AEf 



3fir 



-Pi 



cosfco r f coskor sinfegr 



k r T Hl V (fc r)3 (k ry 



and 



n 3 

r 2 



-Pi 



sin k$r 



sin kgr cos kor 



k r ' V (fc r) 3 
where we have defined ko = ujq/c, 

Pi = 1 - f ? ,qt = 1 - 3 r? 



(12) 

)] 

(13) 
(14) 



and f = (1, 6, ip) being a unit vector along the direction 
joining the two atoms. For a Am = m e — m g = tran- 
sition, 

po = sin 2 0, g = 1 - 3cos 2 6> (15) 
while for a Am = ±1 transition, 

p± = 1(1 + cos 2 9) , q± = 1(3 cos 2 6-1). (16) 

At short distance kor < 1, we obtain that ~ 2r for 
the superradiant state \<p^) = 1 10; 0) and ~ for the 
subradiant state \<f\) = 1 00; 0) . 



denote by V e is obtained from (12) by averaging upon the 
random orientations of the pairs of atoms with respect 
to k. Since, according to (15) and (16), (%) = and 
(pi) = 2/3, we obtain for the average potential V e 



eV e {r) = (AEf) 



fir cos kor 



2 k r 

and the average inverse lifetimes of Dicke states are 



(17) 



(r-> = r l + e 



sin kgr 
k n r 



(18) 



which retain the same features as (13) for kor -C 1. 

Let us now characterize the interaction potential V e . 
Whereas for a single pair of atoms, the potential (12) 
is anisotropic and decays at short distance like l/r 3 , a 
behavior that originates from the transverse part of the 
photon propagator, we obtain that on average over angu- 
lar configurations, the potential (17) between two atoms 
in a Dicke state |Z0) in vacuum becomes isotropic and 
decays like l/r. This behavior coincides with the one ob- 
tained by considering the interaction of two-level atoms 
with a scalar wave. This could have been anticipated 
since in that case the transverse contribution qi to the 
photon propagator averages to 0. A related behavior 
for the orientation average interaction potential has been 
also obtained for the case of an intense radiation field [9] 
and it has recently been investigated in order to study 
effects of a long range and attractive potential between 
atoms in a Bose-Einstein condensate for a far detuned 
light [10]. This latter potential, which is fourth order in 
the coupling to the radiation, corresponds to the inter- 
action energy between two atoms in their ground state 
in the presence of at least one photon. The average po- 
tential V e we have obtained is different from that case: 
it is second order in the coupling to the radiation and 
it corresponds to the interaction energy of Dicke states 
|L0) in vacuum. 



C. Scattering properties 

In order to study the scattering properties of Dicke 
states we introduce the collision operator T{z) = V + 
VG(z)V where V is given by (2) and G(z) = (z - H)- 1 
is the resolvent where the Hamiltonian H is the sum of (1) 
and (2). The matrix element that describes the transition 
amplitude from the initial state \i) = |1 — l;ke), where 
the two atoms are in their ground state in the presence 
of a photon of frequency u> = c|k| and polarization e, to 



the final state |/) 



l-l;k'e' 



B. Average interaction potential 



T={f\T(z = h(u-uo))\i) 



(19) 



For a photon of wave vector k incident on an atomic 
cloud, the potential between two atoms that we shall 



where |k| = |k'|. By using the closure relation we may 
write T as the sum of a superradiant and a subradiant 
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contribution, T = T+ + T~ [12] with 
T± = (flVl^K^lGiz = fcCw-woM^X^W) (20) 

where > arc the Dickc states \L0 > in vacuum. The 
two matrix elements in (20) represent the absorption and 
the emission of a real photon by the pair of atoms. They 
are easily obtained from (2)-(5) and lead to the following 
expressions for the scattering amplitudes 



T+ = Ae5 k - k '> R cos 



k r 



k' r 



CP 



(21) 



and 



T- = Ae^- k '> R sin (*^\ sin (K^) G~ . (22) 



We have defined r = ri r2 , R = (ri + v^)/2 and 

A = ^-d 2 (d -e)(d* ■£'*) 
e il 



(23) 



where the reduced matrix element and the corresponding 
unit vector are 

d=M|M ,d=l( jeme \d\j g m g ). (24) 

The propagators G ± are the expectation values of the 
resolvent in the Dicke states l^), namely G = 
(0 ± |G(M)|</) ± ) where close to resonance 5 = u)—u)o <C wo- 
The propagators result from the sum of an infinite series 
of virtual photon exchanges between the two atoms in 
the pair and are given in terms of (12) and (13) by 



G ± = ( hS - AE ± + ih— 



(25) 



The average propagator is then obtained by averaging 
G^ over the random orientations of the pairs of atoms 
with respect to the wave vector k of the incident photon. 
However, we shall consider in a first stage the effective 
propagator obtained for the case of a scalar wave. This 
accounts to write for the effective propagator the expres- 
sion 



G „ — 



r r e lkor 

6 + i-± 

2 2 k r 



(26) 



where we have used (17) and (18) for the average poten- 
tial and for the average inverse lifetimes. This expres- 
sion constitutes a priori a rough approximation of the 
exact average. We shall calculate later, in section VI, 
the exact expression of the average propagator and show 
that it is rather complicated whereas the approximate 
expression using a scalar wave gives similar qualitative 
results. Therefore, it allows for a better understanding 
of relevant physical quantities such as elastic mean free 
path and group velocity. From now on, we thus use the 



scalar wave approximation in order to provide, in a rather 
simple way, the main features of multiple scattering by 
supcrradiant pairs. 

With the help of (26), the scattering amplitudes are 



T+ = Ae^'^cos i^-) cos (K^) Gt 



(27) 



and 



(±^\sn.(*lL\ G -. (28 ) 



At short distances k$r <C 1, the subradiant amplitude 
T~ becomes negligible as compared to the superradiant 
term T+. Therefore, the potential (17) is attractive and 
decays like 1/r. More precisely, at short distances the 
effective propagator G~ diverges for S/T = l/(2fcor) and 
is purely imaginary for S/T = —l/(2kar). Thus, for 
S/T < l/(2fcor) the imaginary part of the subradiativc 
term (28) is negligible as compared to the imaginary part 
of the superradiative term (27) and for \S\/T < l/(2fcor) 
both the real part and the imaginary parts of (28) are 
negligible as compared to (27). 

We can interpret these results by saying that, at short 
distances (kor C 1), the time evolution of the initial state 



\m) = \em;0) = ^ + ) 



(29) 



corresponds, for times shorter than 1/r, to a periodic 
exchange of a virtual photon between the two atoms at 
the Rabi frequency 



(AE~) - (A£ H 
h 



(30) 



which is much larger than T since, with the help of (17), 



T 

k r 



(31) 



For larger times, the two atoms return to their ground 
state and a real photon (k'e') is emitted. At large dis- 
tances (kgr 3> 1), the Rabi frequency becomes smaller 
than r, so that the excitation energy makes only a few 
oscillations between the two atoms, thus leading to a neg- 
ligible interaction potential. 

We finally notice that the angular distribution of the 
light scattered by two atoms in a superradiant state is 
nearly identical to that of a single atom. This follows 
from the fact that at short distance kor « 1, we can 
neglect higher order multipolar corrections so that the 
corresponding additional phase shift, kor cos??, between 
waves emitted by the two atoms becomes negligible (# 
is the angle between the direction of the emitted photon 
and the axis between the two atoms). 
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IV. COOPERATIVE EFFECTS AND COHERENT 
BACKSCATTERING 

It is interesting to derive the previous results in another 
way that emphasizes the analogy with coherent backscat- 
tcring [1, 2]. To that purpose, we write the scattering 
amplitude T defined previously in (19) as a superposi- 
tion of two "classical", scalar amplitudes, T\ and T 2 [13], 
each of them being a sum of single scattering and double 
scattering contributions, that is 



Ti 



t 



and 



T 2 = 



t 2 Gl 



,i(k-k')T! +t Q o e i(k-rx-k'.r 2 ) 



; i(k-k')-r 3 +t Q Q e i(k-r 2 -k'. ri ) 



Here 



t 



4tt 



k S 



(32) 



(33) 



(34) 



the spatial positions ri and r 2 makes the interference 
term cos(k + k') ■ (ri — r 2 ) vanish in general, with two 
exceptions: 

k + k' ~ : In the direction exactly opposite to the 
direction of incidence, the intensity is twice the classical 
value. This phenomenon is known as coherent backscat- 
tcring. 

ri = r 2 : Closed multiple scattering trajectories which 
are at the origin of the phenomenon of weak localization. 

In (36) the interference term, i.e., the second term in 
the bracket, reaches its maximum value 1 for ri = r 2 so 
that we obtain from (32)- (33) and (27) that T x = T 2 oc 
(1/2)T+, up to a proportionality factor [13]. Thus, the 
total amplitude is given by the superradiant term with 
no subradiant contribution. 



MULTIPLE SCATTERING AND 
COOPERATIVE EFFECTS 



A. Effective self-energy 



is the amplitude of a scalar wave scattered by a single 
atom at the origin and the prefactor i/(l — t 2 G^) where 



Gn — — - 



o ik r 

Airr 



(35) 



accounts for the summation of the series of virtual photon 
exchange between the two scatterers. We can single out 
in the total amplitude T = T\ + T 2 the single scattering 
contribution T s and write the intensity associated to the 
higher order scattering term shown in Figure 1 as 



IT - T„ 



t 2 G 



[1 + cos(k + k') • (n - r 2 )] . 

(36) 




FIG. 1: Schematic representation of the two amplitudes that 
describe double scattering of a scalar wave. The wavy line 
accounts for the exchange of a virtual photon between the two 
atoms. This diagram is analogous to the coherent backscat- 
tering in quantum mesoscopic physics. 

The structure of relation (36) is very reminiscent to 
that of the so-called coherent backscattering intensity 
which occurs in the multiple elastic scattering of light. 
But although they arc analogous, (36) differs from co- 
herent backscattering. In the latter case, averaging over 



We consider now multiple scattering of a photon by 
superradiant pairs built out of atoms separated by a dis- 
tance r and coupled by the attractive interaction poten- 
tial V e . This situation corresponds to a dilute gas that is 
assumed to fulfill 



r < A < rii 1/3 



(37) 



where is the density of pairs and Ao = 2n/ko is the 
atomic transition wavelength. The limiting case (37) cor- 
responds to a situation where the two atoms that form a 
superradiant pair, through exchange of a virtual photon, 
constitute an effective scatterer and cooperative interac- 
tions between otherwise well-separated pairs are negligi- 
ble. Let us stress that we study here a simplified model 
where only pairs of atoms have been taken into account. 
A more realistic model should include higher order terms 
that account for cooperative effects between more than 
two atoms, but we do not consider such higher order 
terms, i.e., including superradiant clusters of three or 
more atoms. The purpose of the current model is to ex- 
amine the contribution of superradiant pairs to the trans- 
port properties of the gas. We use the Edwards model 
[1, 14] to describe the medium as a discrete collection 
of Ni superradiant pairs in a volume Q. Each pair, lo- 
cated at R;, is characterized by its scattering potential 
u(R — R;). Therefore, the disorder potential is given by 



U(R) = ^ «(R - R ( 



(38) 



i=i 



We assume that the scattering potential is short range 
compared to the wavelength, and we approximate it 
by a (conveniently regularized) delta function potential, 
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it(R) = ito<5(R). In the limit of a high density of weakly 
scattering pairs, but with a constant value of t^Uq, it can 
be shown [1] that the correlation function defined by 

B(R - R) = rii J dR"u(R" - R)m(R" - R), (39) 



becomes 



B(R - R) = n lU g(5(R - R). 



(40) 



In other words, in this limit, the Edwards model reduces 
to a Gaussian white noise model characterized by the 
condition (40). 

The Green's function g of a scattered photon is related 
to the free photon Green's function go, i.e., in the absence 
of disorder potential, by the equation [1] 



.9 = go+goUg. 



(41) 



Averaging (41) over disorder and using the properties 
of the Gaussian model discussed above yields the Dyson 
equation 



(g)d = go + go^(g)d 



(42) 



where (■ ■ -)d denotes averaging over the random poten- 
tial. The function S, known as the self-energy, represents 
the sum of all irreducible scattering diagrams. The per- 
tubative expansion of the self-energy in a power series 
controlled by the parameter uiUq is represented in Fig- 
ure 2. 



E = ^ + 

Tji 



x x 



s 2 



x 



£ 4 



FIG. 2: Pertubative expansion of the self-energy in a power 
series in the parameter mUg. Solid lines account for the free 
photon Green's function go. Pairs of dotted lines, connected 
by x, stand for the two-point correlation function B. The 
first term £i, proportional to muo, accounts for independent 
scattering events while the second term £2, proportional to 
nfuQ, describes interference effects between pairs of scatterers. 

For small values of rijUp, the main contribution is ob- 
tained by keeping only the first term Si which describes 
independent scattering events. Therefore, the first con- 
tribution to the self-energy is proportional to the density 
of scatterers and to the average scattering amplitude and 
it is given, for kor <C 1, by 



Si 



where 



-4, 



k Jf ' Je 



12je 

3 2j fl 



(43) 



(44) 



The latter quantity is obtained by averaging A in (27) 
over Zeeman sublevels m that appear in its definition 
given by (23) and (24). 

The additional average, denoted by • • •, is taken over 
distances r up to a maximal value r m which accounts for 
all possible mechanisms that may break the pairs. 

The value of r m can be estimated by comparing the 
kinetic energy if of a superradiant pair to its average 
potential energy . We have K ~ ft 2 / /ir 2 and from 
(17) we obtain that ~ —hT/2kor. Minimizing the 
average energy 



E(r) ~ — 
jir z 



2k r 



with respect to r yields 



k r m = 4—^ 



k r n 



(45) 



(46) 



(47) 



where the speeds v m i n and v m ax have been defined in sec- 
tion II. For typical values, T = 10 7 s _1 and k = 10 7 m _1 
we obtain that kor m ~ 0.05. Thus, we can use the results 
obtained in subsection III-C and consider the superradi- 
ant term only. 

For j g — and j e — 1, A i = 1 and using (26) we 
rewrite (43) as 



Si = 



67m,' 1 



dr 



2k r 



+ i 



(48) 



We stress again that, in our approach, a pair of atoms in a 
superradiant state is considered as a single scatterer and 
the effective medium parameters are derived from Si as it 
will be shown in the next subsections. In contrast to our 
treatment, others [16, 19] consider multiple scattering of 
a real photon by independent atoms and use the second 
term S2, which describes interference effects between the 
scatterers, to calculate corrections to the elastic mean 
free path and to the refractive index of the medium. A 
further comparison between these two points of view is 
given in section VII. 



B. Elastic mean free path 

The elastic mean free path l e is obtained from the 
imaginary part of the self-energy, namely 



— - = — ImSi. 



(49) 



Let us stress that (49) is equivalent, in the case of a dilute 
gas, to the known formula 



1 



n l a e 



(50) 
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where the total cross section, <r e , is obtained for k$r <C 1 
from (27) by means of the optical theorem 



<T, 



2VL 

he 



Im(T+(k = k',e = £'))r 



(51) 



and (•••),« represents an averaging over Zeeman sub-levels. 
The equivalence in this case is proven easily if one uses 
(44) and the usual expression for the inverse lifetime 



r 



d 2 kl 



(52) 



where the reduced matrix clement is defined in (24). 
Therefore, from (48) and (49) we obtain that 



67171, / 5 
-3-/1 k r rn ,- 



where we have defined the function 

,■2 a 



i r u _dx__ 
2w Jo i + {v + 



(53) 



(54) 



The integral is easily carried out analytically and the ex- 
plicit expression is given in Appendix A. It is interesting 
to compare l e to the elastic mean free path Iq that corre- 
sponds to near resonant elastic scattering of a photon by 
a single atom. The latter quantity is obtained by replac- 
ing r by r/2 in (53) (since the inverse lifetime of a single 
atom is half the one related to a superradiant pair) and 
1/x by in (54) (since the inter-atomic distance is taken 
to be infinite for a single atom) and it is given by 



lo(S) = 



k 2 
ft-o 

67m; 



(55) 



In Figure 3 the ratio between these two quantities is plot- 
ted as a function of the reduced detuning S/T from res- 
onance for several values of kor m . 

At resonance, we obtain from (53) that 



le(0) 



k 2 



1 



and hence 



*o(0) 

Ze(0) 



= o^m) 2 < I- 



(56) 



(57) 



Away from resonance, for blue detuning, the elastic mean 
free path l e becomes smaller than ^ in a ratio roughly 
given by l/(/cor m ) 2 . This is a direct consequence of the 
existence of the attractive potential V e . 



C. Group velocity 

Another important physical quantity that character- 
izes multiple scattering of a photon is its group velocity 




FIG. 3: Ratio between the elastic mean free paths lo and l e as 
a function of the reduced detuning S/T for kor m = 0.05, 0.07, 
and 0.1. Away from resonance, for blue detuning, the elastic 
mean free path l e becomes smaller than Iq in a ratio roughly 
given by l/(kor m ) 2 . At resonance, the ratio between the elas- 
tic mean free paths is given by (57). 



v g given in terms of the refraction index 77 by the usual 
relation 



c drj 
— = ri + uj—. 



(58) 



The refraction index for a dilute medium is 

r) = (1 +n 4 Rea) 1/2 (59) 

where the dynamic atomic polarizability a is propor- 
tional to the self-energy 



a = --(-) 

rii \uis 



Thus, we obtain that 



-1 1/2 



RcS] 



Substituting (61) into (58) yields 



- = - 1 - — — ReSi 
v g 77 \ Zlu aw 



(60) 



(61) 



(62) 



From the self-energy (48), we notice that 77 ~ 1 for all val- 
ues of the detuning S /T and in a large range of densities 
71, so that 



Vg(S) 



Hi ( S 

1 72 k a r rn , - 

n r _ \ I 



where we have defined the characteristic density 

67T CJQ 



and the function 



1 



dx- 



(63) 



(64) 



(65) 
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The integration is easily performed and the explicit ex- 
pression is given in Appendix A. By replacing T by T/2 
in (63) and 1/x by in (65), we obtain the group velocity 
Vq of light interacting with independent two level-atoms 



v (S) 



Mi + (¥) 2 ) 



(66) 



For the following typical values, T 



10 7 i 



10 m and 



10 



10, 



, we obtain that rii/n c ~ 10 



Figure 4 displays the group velocities v g and vq plotted 
as a function of the reduced detuning S/T for rii/n c = 10° 
and fco7' m = 0.1. 




FIG. 4: Group velocities v a (solid line) and vo (dotted line) 
as a function of the reduced detuning 5/T for ni/n c = 10 
and kor m = 0.1. The group velocity vo diverges at two sym- 
metric values of order unity of the reduced detuning and it 
takes negative values in between. The group velocity v g , near 
resonance, remains finite and positive. 

v g appears to diverge at quite a large and negative 
value of the detuning S/T ~ — l/(2fc r m ). But near res- 
onance it is well behaved, meaning that it remains finite 
and positive. At resonance, according to (63), the group 
velocity is 



c , . . rii uj , , 2 

^ = 1 + 4 ^T- (fcorm) - 



(67) 



This expression of v g differs substantially from the one 
obtained for vq- For densities n, > n c , the group velocity 
vq diverges at two symmetric values of order unity of the 
detuning and it takes negative values in between (i.e., 
also at resonance), as it can be seen in Figure 4. This 
problem has been recognized a long time ago [15] and an 
energy velocity has been defined which describes energy 
transport through a diffusive medium [16, 17]. However, 
the diffusion coefficient, which will be discussed in the 
next subsection, is derived from the group velocity and 
not from the energy velocity [1]. Moreover, a closed ex- 
pression for the energy velocity, ve, has been obtained 
only for the case of a resonant Mie scattering [18]. The 



expression is similar to (67) and is given by 
c niUJo 

= 1 + 97T — 

VE 



kl r 



(68) 



It is then interesting to notice that the inclusion of coop- 
erative effects even at the lowest order, i.e., taking into 
account superradiant pairs, allows to obtain a group ve- 
locity which is well behaved at resonance, unlike the case 
of resonant scattering by independent atoms. 



D. Diffusion coefficient and transport time 

Diffusive transport of photons through a gas is charac- 
terized by the photon diffusion coefficient 



D(S) = -v g (S)l e (S) 



(69) 



that combines the elastic mean free path and the group 
velocity, both derived from the complex valued self- 
energy (48). The diffusion coefficient D is of great im- 
portance since it enters into expressions of various mea- 
sured physical quantities, such as the transmission and 
the reflection coefficients of a disordered medium [1]. In 
addition to these average quantities, an incident pulse 
that probes a nearly static configuration of scatterers 
may provide an instantaneous picture of the medium 
that displays a random distribution of bright and dark 
spots. This snapshot, known as a speckle pattern, can 
be characterized by the angular-correlation function and 
the time-correlation function of the light intensity (diffus- 
ing wave spectroscopy). In the first case, the correlation 
function of the transmission coefficient between two dis- 
tinct directions of the transmitted wave is measured. In 
the second case, the intensity of the transmitted wave 
is measured at different times, so that the motion of 
the scatterers must be taken into account. As pointed 
before, in both cases the diffusion coefficient plays an 
important role, as it enters in the relevant expressions. 
Moreover, the critical behavior of transport close to the 
Anderson localization transition at strong disorder is also 
obtained from the scaling form of D. Its expression, de- 
duced from (53) and (63), depends on the range r m and 
on the detuning S/T. Since the group velocity and the 
clastic mean free path are significantly modified for su- 
perradiant states, we thus expect the diffusion coefficient 
to be different from its value obtained for independent 
atoms. 

We define the transport time by 



Ttr{5) 



US) 



(70) 



At resonance and for ri; > n c , it can be rewritten with 
the help of (56) and (67) as 

1 

2F 



rtr(0) 



(71) 
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in accordance with our assumption of supcrradiant 
states. Near resonance, the transport time depends 
weakly on the detuning. But, away from it, r tr depends 
on the detuning and thus on frequency, as it can seen 
from Figure 5 where the inverse of the transport time 
T t ^ 1 /r is plotted as a function of the reduced detuning 
5/T for m = 10 10 cm- 3 , T = 10 7 s _1 and fc = lOW 1 
for several values of kor m . 




sir 



FIG. 5: Inverse of the transport time r t ~ 1 /r as a function of 
the reduced detuning 5/T for m = 10 10 cm -3 , F = 10 7 s~ 1 and 
ko — \Q 7 m~ 1 for kor m = 0.05, 0.07, and 0.1. Near resonance, 
the transport time depends weakly on the detuning. But, away 
from it, Ttr depends on the detuning and thus on frequency. 



VI. AVERAGE SELF-ENERGY 



In Figure 6 the ratio between l Q given by (55) and l' e 
is plotted as a function of the reduced detuning 5/T for 
several values of kor m . 




-100 -80 -60 -40 -20 

8/r 



FIG. 6: Ratio between the elastic mean free paths lo and l' e as 
a function of the reduced detuning 5/T for kor m = 0.05, 0.07, 
and 0.1. At resonance l' e is larger than lo, but away from 
resonance it becomes smaller. 

As in the effective approach, at resonance l' e is found 
to be larger than Z , but away from resonance it becomes 
smaller. 

In Figure 7 the group velocity v' g is plotted as a func- 
tion of the reduced detuning 5/T for ni/n c = 10 5 and 
k r m =0.1. 



So far, we have used the effective approach introduced 
in subsection III-C, where we have considered the case of 
a scalar wave being scattered by a pair of two-level atoms. 
In this simple approach, the propagator of a scalar wave 
(26) has been calculated and the self-energy (43) has been 
obtained by averaging (26) over the distance between the 
two atoms in a pair. This effective approach leads to 
simple expressions for the elastic mean free path (53) and 
the group velocity (63) of the wave. In this section we 
calculate these quantities for a given Am transition and 
kor <C 1, while taking into account the vectorial nature 
of the wave. To this purpose, we average the propagator 
(25) over the random orientations of the pairs of atoms 
(with respect to the wave vector of the incident photon) 
as well as over the distance between the two atoms in a 
pair. Therefore, the average self-energy is now given by 

= f hTG+(r)dr (72) 

fc 4nr m J 

where the averaging is over the inter-atomic axis r (both 
over magnitude and orientations). The evaluation of (72) 
for a Am = transition is rather cumbersome and it is 
presented in Appendix B. By following the procedure de- 
scribed in the previous section, we obtain the correspond- 
ing elastic mean free path l' e and the group velocity v' g . 



1.0086 
1.0085 
1 .0084 
1.0083 
1 .0082 
1 .0081 
1.008 
1.0079 
1.0078 
1 .0077 



8/r 

FIG. 7: Group velocity v' g as a function of the reduced detun- 
ing 5/T for m/ric = 10 5 and kor m — 0.1. Around resonance, 
the group velocity v g is finite and positive and it is close to c. 

Around resonance, the group velocity v' g is finite and 
positive, as in the scalar case, but much larger as com- 
pared to (63) and it is close to c. Thus, we may conclude 
that in both approaches the superradiant effect leads to a 
finite and positive group velocity, unlike the one obtained 
for light interaction with independent atoms. However, 
the group velocity of a scalar wave is much smaller com- 
pared to the one of a photon. 
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VII. DISCUSSION 



APPENDIX A 



In this section we compare our analysis to other ap- 
proaches [16, 19] where resonant multiple scattering of 
light has been considered. There, using a multiple scat- 
tering expansion for the calculation of the self-energy up 
to second order in run 2 ,, a correction to the elastic mean 
free path and to the refractive index has been obtained. 
In the latter approach, no distinction has been made be- 
tween the external photon that performs multiple scat- 
tering on all atoms and virtual photons exchanged be- 
tween two atoms in a superradiant state leading to the 
average interaction potential V e . This distinction needs 
to be made for dilute enough atomic gases since in that 
case the average distance n i 1//3 between atoms is large. 
Moreover, in this case, the dipolc-dipolc interaction in- 
duced by the external photon depends on the detuning, 
a situation that corresponds to the case of intense radia- 
tion presented in [9] but not to the current experiments 
made on cold atomic clouds [20]. 



In this Appendix, we establish expressions (53) and 
(63) for the elastic mean free path and the group veloc- 
ity. At resonance, simple expressions for the elastic mean 
free path (56) and the group velocity (67) are obtained 
by a pertubative expansion with respect to the small pa- 
rameter kor m . 



1. Elastic mean free path 

The elastic mean free path is given by (53) in terms of 
the function /i defined in (54). The integral in (54) is 
easily carried out analytically and it leads to 



1 



67m, 1 



l e (5) k 2 ad \2a 
where 



b A 1 + a -^B 1+Cl 



2a 



a = l + {S/T) 2 b = S/T 



(Al) 



(A2) 



VIII. CONCLUSIONS 

We have considered multiple scattering of a photon 
by pairs of atoms that are in a superradiant state. On 
average over disorder configurations, an attractive inter- 
action potential builds up between close enough atoms, 
that decays like 1/r. The contribution of superradiant 
pairs, resulting from this potential, to scattering prop- 
erties is significantly different from that of independent 
atoms. This shows up in the behaviors of the group veloc- 
ity, the elastic mean free path and the diffusion coefficient 
which arc different from their values obtained for inde- 
pendent atoms. We have considered the case of a scalar 
wave and have shown that it allows to define an effective 
long range and attractive potential for pairs of atoms in 
a superradiant state. Then, we have studied the case of a 
vector wave and have shown that the results obtained in 
the scalar case remain qualitatively valid. We have con- 
sidered a simplified model where only pairs of atoms have 
been taken into account. A more realistic model should 
include higher order terms that account for cooperative 
effects between more than two atoms [21]. The purpose 
of the current model is to show that already for a dilute 
gas in the weak disorder limit, cooperative effects modify 
significantly the transport properties of light. 



and 



A l = In 



±x 2 



d -\- bx m -\- 



7T _, 1 

Bi = - - tan (6 + -x m ) 



k r m < 1. 



(A3) 



(A4) 



(A5) 



At resonance (5 = 0) we have a = 1, b = and by 
expanding (A4) with respect to fco r m we obtain 



Thus, 



2 / 4 



k 2 1 

ie(0) - 



as given in (56). 



8TTrii (kor m y 



2. Group velocity 



(A6) 



(A7) 
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The group velocity is given by (63) in terms of the 
function defined in (65). The integral in (65) is easily 
carried out analytically and it yields 



v g (S) 



= 1 - 



Fi 



n c a 2 C\ 



(A8) 



where 
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^ = 6(~^ 1+ ^i + (|-f)fl 1 - fc fli + (l-f)C 1 

(A9) 



(1 ~\- bx m -\- 



and 



i 

2 



(A10) 



(All) 



At resonance (S = 0) we have a = 1, b = and by ex- 
panding (A10) and (A4) with respect to k r m we obtain 



and 



Thus, 



2/4 
A i ~ — — - 1 

2-m V^m 



2 / 4 



w fl (0) 



= 1 + - — (k a r„ 
6 n c 



(A12) 



(A13) 



(A14) 



as given in (67). 



APPENDIX B 



The aim of this Appendix is to calculate the average 
self-energy (72) for a Am = transition in the case where 
kgr <C 1. First, we average the supcrradiative propaga- 
tor (25) over the orientation of the inter-atomic axis and 
obtain analytical expressions for its real and imaginary 
parts. Then, by averaging over the inter-atomic distance 
up to r m , wc obtain the average self-energy (72). 

For a Am = transition and k$r <C 1, the superra- 
diative propagator (25) may be written with the help of 
(12) and (13) as 



hTG + = 



6 3 /3cos 2 6» - 1 , |(1 + cos 2 i 
f + " 



(fc r) 2 



fc r 



+ i 



(Bl) 

where the inter-atomic axis is r = (r, 9, ip). An averaging 
over the orientations 



nr(G + ) = — [ hTG+dcos 
4vr J 

yields for the imaginary part, 

P + Q 



hTlm(G + ) 



(B2) 



(B3) 



and for the real part, 

hTRe(G+) = W-P + W + ( 
where we have defined 



P 



1 



-hx 



Q 



1 + 2/3cos(f) +0- 
8A 2 /3cos(fp \l - 2/3 cos(%) + /3- 

-i 1-/3 2 



4A 2 /3sin(f) V2 



tan 



and 



W± = -v^2(cos7 =F !)• 
The auxiliary parameters are given by 



0=& 7 = 008-^-5^ 



where 



.4, 



16(fc r) 2 V(M 



(B4) 

(B5) 
(B6) 

(B7) 
(B8) 

(B9) 



B 2 = 
and 



C 2 = 1 + - 





3 




28 


3 




(fcor) 2 




r 


2fc r 




1 


'25 


3 


a 


1 




4 


r 4 


2fc r 




(fcor 





1 1 



(BIO) 
(Bll) 



Finally, we average (B3) and (B4) over the inter-atomic 
distance up to r m 

hTlm(G+) = -— / dr^f- (B12) 

r m Jq P 

and 



hTRc(G^ 



' m JO 



dr(W-P + W+Q). (B13) 



The integrals can be evaluated numerically and give the 
average self-energy (72) since 



1 



47rr I7 



frTG + (r)dr = hT(G+). (B14) 
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